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ABSTRACT 

An  optimal  integral  control  design  program  was  applied 
to  a  linearized  state  variable  model  of  the  FFG-7  ship  class 
gas  turbine  and  Controllable  Reversible  Pitch  (CRP)  propeller 
main  propulsion  system.   Various  combinations  of  output 
parameters  were  investigated  in  an  attempt  to  produce  a 
feasible  control  design.   Only  one  acceptable  design  was 
achieved  which  did  not  violate  any  physical  constraints. 
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I.   INTRODUCTION 

A.  OBJECTIVE 

Advanced  technology  propulsion  systems  require  a  control 
system  design  that  is  able  to  achieve  the  desired  performance 
at  the  greatest  possible  efficiency.   One  particular  scheme 
that  shows  great  promise  in  this  regard  is  the  optimum  integral 
control  design  program,  CONSYN,  developed  in  Ref.  1. 

This  study  was  undertaken  to  examine  the  validity  and 
effectiveness  of  applying  CONSYN  to  the  FFG-7  ship  class  gas 
turbine  propulsion  system  model  obtained  from  Ref.  2. 

B.  COMPUTER  FACILITIES  AND  SOFTWARE  UTILIZED 

The  above  objective  was  implemented  by  simulation  on  the 
IBM  360/67  computer  of  the  W.  R.  Church  Computer  Center  at 
the  Naval  Postgraduate  School.   The  following  programming 
software  was  required: 

- 

Various  routines  in  the  International  Mathematical  and 
Statistical  Library  (IMSL) 

The  IBM  developed  Continuous  System  Modeling  Program 
(CSMP) 
*  The  Control  Program  for  Engineering  Synthesis  with 
Constrained  Function  Minimization  (COPES/CONMIN)  de- 
veloped by  Garret  N.  Vanderplaats  of  the  NASA  Ames 
Research  Center 
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The  routines  developed  for  CONSYN  in  Ref.  1  as  modified 
by  Ref.  3  and  this  study,  are  listed  in  Appendix  A. 


II.   PROPULSION  PLANT  MODEL 

The  propulsion  system  for  the  OLIVER  HAZARD  PERRY  (FFG-7) 
class  patrol  frigates  consists  of  two  General  Electric 
Corporation  LM2  500  Marine  Gas  Turbines,  coupled  through  a 
reduction  gear  to  a  controllable  reversible  pitch  (CRP) 
propeller  manufactured  by  the  Bird-Johnson  Company.   Since 
the  gas  turbine  is  a  highly  non-linear  device,  the  linearized 
state  variable  model  developed  in  Ref.  2  is  valid  only  for 
small  changes  from  a  specified  steady  state  operating  point. 
This  model  was  obtained  from  a  coupled  ship  and  propulsion 
plant  dynamic  non-linear  model  consisting  of  the  following 
features : 

*  A  non-linear  model  of  the  LM2  50  0  engine  developed  by 
General  Electric. 

A  non-linear  model  of  the  ship,  propeller  and  propul- 
sion train  provided  by  the  Bath  Iron  Works  Company. 
Secondary  flow  losses  considered. 

Bleed  flow  for  extraction  and  turbine  cooling  considered. 
Inlet  and  exhaust  losses  incorporated. 
Pressure  and  temperature  dynamics  neglected. 
Nested  loop  balancing  utilized  to  obtain  various  steady 
state  operating  points. 

In  matrix  equation  form,  this  linearized  state  variable 
model  of  the  various  engine  and  ship  parameters  is  expressed 
as : 
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X  =  AX  +  BU 
Y  =  CX  +  DU 
where 

X  is  the  system  state  vector,  and  consists  of  three 
components : 

NG  =  Gas  Generator  Speed  (rpm) 
NPT  =  Power  Turbine  Speed  (rpm) 
V  =  Ship  Speed  (knots) 


X-,  =  NG  -  NG_ 

1  ss 

x_  =  NPT  -  NPT 

2  ss 


x   =  V  -  V 
3         ss 


U  is  the  system  input  vector,  and  consists  of  four 
components : 


u,  =  WF  -  WF 
1  ss 


u,  =  3  -  300 

2.  SS 

u.,  =  WB  -  WB 

3  ss 

u.  =  P/D  -  P/D 

4  '     ss 


WF  =  Fuel  Flow  (lbm/hr) 
3  =  Stator  Guide  Vane  Angle  (deg) 
WB  =  Compressor  Bleed  Flow  (lbm/s) 
P/D  =  Propeller  Pitch/Diam  Ratio  ('*') 
Y  is  the  system  output  vector,  and  consists  of  six 
components : 

SMC  =  Compressor  Surge  Margin  (%) 

T4  =  Combustor  Outlet  Temp  (°R) 

QPT  =  Power  Turbine  Torque  (ft-lbf) 

P3  =  Compressor  Total  Disch  Press  (psia) 

HPPT  =  Power  Turbine  Power  (hp) 


y,  =  SMC  =  SMC 

y0  =  T4  -  T4 
1  2  ss 

y3  =  QPT  -  QPT 

y.  =  P3  -  P3 
J4  ss 


ss 


ss 


y5  =  HPPT  -  HPPT 

y,  =  T5  -  T5 
2  6  ss 


ss 


T5  =  Power  Turbine  Inlet  Temo  (  R) 


and  the  ss  subscripts  indicate  values  at  a  specific  steady  state 
operating  point,  the  same  point  used  for  determination  of  the 
constant  A,  B,  C,  D  matrix  elements. 

This  study  considered  single  engine  up-power  transients  of 
70%  to  7  5%,  and  8  0%  to  8  5%  of  maximum  gas  generator  speed.  In 
both  cases,  the  steady  state  operating  point  selected  was  that 
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for  the  higher  power,  or  final  condition.   The  necessary 
state  variable  data  to  accomplish  this  study  was  obtained 
from  computer  printouts  generated  in  conjunction  with  Ref.  2, 
and  is  contained  in  Appendix  B.   Physical  constraints  on 
inputs  and  input  rates  are  also  contained  in  Appendix  3. 
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III.   CONTROL  SYSTEM  DESIGN  WITH  CONSYN 

A.   BASIC  PHILOSOPHY 

The  CONSYN  (Control  Synthesis)  design  program  presented 
in  Ref.  1  is  a  further  refinement  of  linear  regulator  theory 
developed  in  Ref.  4.   It  is  based  on  the  initial  assumption 
that  the  plant  dynamics  can  be  represented  by  the  linearized 
state  variable  matrix  equations 

X  =  AX  +  BU 
Y  =  CX  4-  DU 

where,  as  previously  discussed,  all  states  are  perturbed 
from  their  steady  state  operating  points.   This  is  illus- 
trated in  block  diagram  form  in  Fig.  la.   This  original  plant 
is  then  augmented  with  integrators  to  achieve  open  loop  inte- 
gral control  as  shown  in  Fig.  lb.   Application  of  linear 
quadratic  regulator  theory  to  this  augmented  system  leads 
to  the  minimization  of  a  performance  or  cost  index,  J, 
defined  by  the  matrix  equation 

fcf 
J  =  /    (YTQY  +  UTRU)  dt 
o 

where  Q  is  a  constant,  diagonal,  positive  semi-definite 
matrix  used  to  weight  output  deviation  from  final  values,  and 
R  is  a  constant,  diagonal,  positive  definite  matrix  used  to 
weight  input  rate  deviations  from  final  values . 
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The  linear  regulator  problem  solution  yields  an  optimal 
state  feedback  regulator  which  follows  the  matrix  equation 
control  law 

U  =  Z  -  GLX  -  G2U 

and  is  illustrated  in  Fig.  lc  ,  where  Z  is  the  demand  vector. 
Finally,  in  order  to  achieve  integral  control  commands  gen- 
erated by  the  error  between  demand  and  output,  the  system 
shown  in  Fig.  lc  is  converted  to  that  shown  in  Fig.  Id.   This 
final  system  configuration  has  the  same  time  response  as  that 
shown  in  Fig.  lc  but  has  the  added  advantage  of  assuring  zero 
steady  state  error  even  if  system  degradation  alters  the 
elements  of  coefficient  matrices. 

Since  that  choice  of  Q  and  R  matrix  elements  for  the 
solution  of  the  linear  regulator  problem  is  arbitrary, 
COPES/CONMIN  logic  is  then  utilized  to  find  the  optimum  Q  and 
R  which  minimizes  the  cost  index  J  without  violating  any  of  the 
imposed  physical  constraints. 

In  summary  then,  the  CONSYN  design  program  applies  a 
linear  regulator  design  form  to  the  given  output  control 
problem  using  the  augmented  state  vector 

Block  diagram  methods  along  with  suitably  partitioned  matrices 
are  utilized  to  arrive  at  the  desired  structure  of  the  optimal 
control  system. 
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B.   PROGRAM  FLOW  LOGIC 

The  CONSYN  design  program  operates  in  three  distinct 
segments:   Initial  Design,  Optimization,  and  Final  Design 
Output.   A  summary  of  the  significant  logic  and  operations 
within  each  of  these  segments  follows. 

1 .   Initial  Design 

a.  Input  the  following: 

(1)  A,  B,  C,  D  matrices 

(2)  Initial  guess  for  appropriate  Q  and  R  matrices 

(3)  Initial  condition  X  and  U  vectors,  which  are 
then  combined  to  form  the  augmented  state 
vector   X*  =  [-jy-] 

(4)  System  constraints 

b.  Calculate  system  open  loop  characteristics 

c.  Calculate  P,  the  solution  to  the  steady  state 
matrix  Riccati  equation 

A*  P  +  PA*  -  PB*  R"1  3^  P  +  cl    Q  C*  =  0 

where  A*  =  [-§]■§-]    B*  =  [-§-]    c*  =  [C  j  Dj 

i  i 

d.  Calculate  G  the  optimal  state  feedback  gain  matrix 

G  =  R_1  3*  P  =  [G]_  !g2] 


The  detailed  development  of  this  Riccati  equation  and 
the  following  G,  H,  and  L  gain  relationships  in  the  partitioned 
matrix  form  is  found  in  Refs.  1  and  4. 
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e.  Calculate  H  and  L  the  integral  control  gain  matrices 

[l|h]  -  o  c-J-fj-]  _1 

i 

where  E  =  [i  iO]  and  is  required  so  that  the  number 
i 

of  outputs  fed  back  for  comparison  with  the  demand 
vector  is  equal  to  the  number  of  inputs. 

f.  Simulate  time  response  of  system. 

g.  Print  results  of  all  of  the  above. 

2 .   Optimization 

Vary  Q  and  R  matrix  elements  under  the  direction  of 
COPES/CONMIN  logic  to  achieve  the  smallest  cost  function 
value  without  violating  any  constraints. 
3  .   Final  Design  Output 

Print  complete  results  of  system  characteristics 
associated  with  the  optimum  Q  and  R  matrix  elements. 
Since  the  CONSYN  program  does  not  contain  a  plotting  routine 
for  system  time  response  simulation,  a  separate  program  must  be 
utilized.  For  this  study,  the  IBM  CSMP  in  conjunction  with  a 
VERSATEC  plotter  was  employed.   A  derivation  of  the  equations 
required  for  this  simulation  is  contained  in  Appendix  C. 

C.   MODIFICATIONS  TO  THE  ORIGINAL  PROGRAM 

As  discussed  in  Ref.  3,  the  following  two  significant  modi- 
fications were  made  to  the  original  CONSYN  program: 

The  Kleinman  technique  for  the  solution  of  the  matrix 

Riccati  equation  in  subroutine  OPGAIN  was  replaced  by  an 

eigenvalue  method. 
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The  fixed  step  integration  solution  in  subroutine  PEAK 
for  system  simulation  was  replaced  by  a  variable  step 
integration  routine. 

In  using  CONSYN  with  these  changes  for  this  study,  diffi- 
culty was  encountered  with  the  variable  step  integration  routine, 
Numerical  underflows,  overflows,  and  step  size  decimation  were 
encountered  to  varying  degrees  in  all  attempts  to  use  it.   Con- 
sequently, the  fixed  step  integration  method  was  returned  to 
subroutine  PEAK,  along  with  a  logic  modification  to  determine 

and  print  out  the  number  of  terms  required  for  the  matrix 

Ft 
series  convergence  of  PHI  -  e   .   Here,  the  matrix  F  is 

defined  as 

F  =  A*  -  B*  G 

where  the  matrices  A* ,  B*  and  G  are  as  previously  defined. 
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Figure  1.   Integral  Control  Synthesis 
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Fig.  lc .   Closed  Loop  Optimal  Integral  Control 
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Figure  1.   (Cont'd; 


19 


IV.   RESULTS 

A.   GENERAL 

The  overall  results  were  disappointing.   Of  all  the 
basic  design  forms  and  variations  evaluated,  only  one  was 
found  which  did  not  violate  one  or  more  physical  constraints. 
In  addition,  the  study  was  hampered  by  the  relatively  long 
computer  execution  times  required  for  the  optimization  of 
a  single  proposed  design  (approximately  20  minutes) ,  and  the 
excessively  long  computer  turn-around  times  (approximately 
10  hours) .   In  all,  at  least  150  difference  designs  were 
attempted,  with  a  total  computer  execution  time  expenditure 
of  over  50  hours.   A  summary  of  the  design  forms  and  varia- 
tions considered  is  given  in  Tables  I  and  II. 

Initial  efforts  to  find  an  acceptable  design  were  con- 
centrated on  design  form  numbers  2  and  4  in  Table  I.   These 
forms  were  selected  because  they  appeared  to  be  the  most 
analogous  to  the  design  evaluated  in  Ref .  5  for  the  F4  01 
turbofan  jet  aircraft  engine  with  respect  to  those  outputs 
designated  for  feedback.    All  proposed  designs  of  these 
forms  consisted  of  transients  between  70%  and  7  5%  of  maximum 
gas  generator  speed,  variable  bleed,  and  no  state  variable 
scaling.   All  variations  considered  resulted  in  violation 


See  Appendix  D  for  a  more  detailed  discussion. 
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of  one  or  more  physical  constraints.   Examples  most  often 
noted  were : 

Maximum  propeller  pitch  exceeded 

*  Maximum  propeller  pitch  rate  exceeded,  either  in  the 
positive  or  negative  direction 

Negative  bleed 

Design  form  number  4  was  then  exercised  at  various 
transients  between  30%  and  85%  of  maximum  gas  generator 
speed  with  similar  types  of  physical  constraint  violations. 

In  an  effort  to  determine  why  negative  bleed  was  some- 
times occurring,  G,  J.  Michael,  one  of  the  authors  of  Ref s . 
4  and  5  was  contacted.   The  following  information  and  guidance 
was  obtained  from  this  communication: 

Although  variable  bleed  was  used  as  a  control  in  Ref.  4, 
it  also  sometimes  assumed  negative  values  due  to  mathe- 
matical anomalies  in  the  engine  model  that  could  not  be 
eliminated.   As  a  result,  bleed  was  not  used  as  a 
control  input  in  Ref.  5. 

Although  the  primary  motivation  for  state  variable 
scaling  in  Ref s .  4  and  5  was  to  avoid  publishing 
proprietary  information  in  the  open  literature,  doing 
so  in  the  case  of  this  study  might  eliminate  additional 
computational  anomalities. 
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u  = 


Y  = 


TABLE  I 
DESIGN  FORMS  EVALUATED 


Basic  Form  I.D.  Number 


2 

4 

5 

6 

7 

3 

Ul 

Ul 

Ul 

ul 

Ul 

Ul 

U2 

U2 

U2 

U2 

U2 

U2 

U3 

U3 

U4 

U4 

U4 

U4 

U4 

U4 

U2 

U2 

U2 

U2 

U2 

Ul 

U4 

x3 

X3 

y3 

y5 

U4 

y3 

y3 

y6 

y6 

^e 

y6 

X3 

y6 

y3 
yl 

X3 

yl 

x3 
*1 

*1 

*1 

y± 

y2 

y2 

y2 

y2 

y2 

y2 

y3 

y4 

y4 

y4 

y4 

y4 

y5 

Y5 

Y5 

y5 

Y5 

y3 

y6 

NOTE:   The  first  n  elements  of  the  Y  vector,  where  n  is  the 
number  of  inputs,  are  those  that  are  to  be  fed  back 
for  comparison  with  the  demand  vector. 


22 


TABLE  II 
DESIGN  FORM  VARIATIONS  EVALUATED 

1.  Initial  to  final  Power  Level,  %  of  Maximum  Gas  Generator 
Speed : 

a.  70  -  75% 

b.  80  -  85% 

2.  Bleed: 

a.  Variable 

b.  Constant 

3.  Settling  time  -  20  to  200  seconds. 

4.  Ship  speed  constraint  with  respect  to  settling  time: 

a.  Within  +  2%  of  steady  state  value 

b.  Not  constrained 

5.  State  variable  scaling  or  normalization  with  respect  to 
the  difference  between  initial  and  final  steady  state 
values : 

a.  Scaled 

b.  Unsealed 

6.  Elements  of  the  Q  and  R  performance  index  weighting 
matrices  allowed  to  be  changed  by  the  optimization 
routine : 

a.  All  Q  and  R  elements 

b.  r,  (the  weight  of  fuel  flow  rate  error)  fixed  at 
1.0,  all  other  Q  and  R  elements  allowed  to  change, 
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Design  form  numbers  5,  6,  7  and  8  were  then  exercised 
at  various  transients  between  80%  and  85%  of  maximum  gas 
generator  speed,  constant  bleed,  and  with  state  variable 
scaling.    The  following  results  were  noted: 

An  acceptable  design  which  did  not  violate  any  of  the 
input  or  input  rate  constraints  was  achieved  with 
design  form  number  6,  with  the  performance  index 
weight  on  fuel  flow  rate  error  fixed  at  1.0. 

There  was  no  significant  difference  in  system  response 
obtained  from  design  form  numbers  5,  6  and  7. 

For  those  designs  which  did  not  impose  a  specification 
on  ship  speed  with  respect  to  requested  settling  time, 
the  optimization  routine  reduced  the  degree  of  initially 
violated  constraints  during  the  first  several  intera- 
tions .   The  remaining  iterations  produced  no  further 
changes  in  the  degree  of  constraint  violation  or  in 
the  value  of  the  objective  function. 

For  those  designs  which  did  impose  a  specification  on 
ship  speed  with  respect  to  requested  settling  time, 
the  degree  of  initially  violated  constraints  was 
reduced  at  the  expense  of  violating  another  constraint 
that  was  not  initially  violated. 


See  Appendix  E  for  the  derivation  of  the  scaled  state 
variable  data  with  constant  bleed. 
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The  number  of  terms  required  for  the  series  convergence 

Ft  . 
of  PHI  =  e    in  subroutine  PEAK  was  much  less  than 

before  state  variable  scaling  was  utilized;  usually 

fewer  than  ten  terms  were  needed  for  convergence. 

Deviations  from  actual  values  were  noted  for  the 
linearized  model  calculations  for  outputs  at  the  initial 
condition.   The  worst  occurred  for  power  turbine  power 
which  assumed  a  negative  value. 

B.   SPECIFIC 

A  summary  of  the  design  data  associated  with  the  one 
acceptable  design  achieved  is  presented  in  Table  III.   Plots 
of  system  time  response  simulation  for  this  design  are  shown 
in  Figures  2a  through  2f.   Plots  of  state  open  loop  response 
to  a  step  input  corresponding  to  the  final  steady  state  input 
values  are  shown  in  Figures  3a  and  3b.   It  is  noted  that  with 
the  optimum  design,  all  engine  parameters  reached  their  final 
steady  state  values  within  6  seconds,  and  the  ship  speed 
achieved  99%  of  its  final  steady  state  value  after  100  seconds 

Summaries  of  performance  data  for  various  selected  designs 
are  presented  in  Table  IV,  with  the  degree  of  constraint 
violation  indicated. 

A  comparison  of  actual  output  values  with  those  predicted 
by  the  linearized  model  at  the  initial  condition  is  given  in 
Table  V. 
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TABLE  III 
OPTIMAL  CONTROLLER  DATA,  DESIGN  6-2F 
Performance  Index  Weighting  Factors 

Output  Error 


(u2)  (y3: 


66 


.37 


(y6)   (x3)   (yL) 


26   7x10 


=7 


6x10 


°*6 

(y2) 

.56 


^7 

(y4) 

.72 


(y5) 


Input  Rate  Error 


(ux) 


023  I  1.0 


"2 

:u2j 

.32 


(u4) 

.23 


The  H  Matrix 


.3032 

.5066 

.9441 

1.533 

.03452 

-.02609 

-.1797 

.9716 

-.3782 

The  L  Matrix 


-.3585 

-.02830 

-.0006817 

.1184 

-.03476 

-.001542 

.8133 

-.3113 

-.2850 

NOTE:   The  complete  computer  output  of  the  optimization  run  for 
this  design  is  contained  in  Appendix  F . 
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TABLE  IV 
PERFORMANCE  DATA  SUMMARY  FOR  VARIOUS  DESIGNS 

Design  Code 


Constraint   Limit    2-lA   4-lA  4-1B  4-2B  6-2C   6-2D   6-2E   6-2F   3-2?:   3-2F 

1640   1640  1350  3360  3370   4170   3860   3360   3950   3910 

1320   1320  1320  2360  2360   2360   2360   2360   2360   2260 

33.3   33.3  33.3  23.4  23.4   23.4   23.4   23.4   23.4   23.4 

29.0   27.5  29.4  13.0  13.3   15.0   13.2   13.2   11.9   13.1 

1.61   2.54  2.03  10.4  constant  =  .02 

JTTT9)  ^TTT|>  .02  (£  ~ 22) constant    =    .02 

^2~To~5)   (2TT5)  (2TT0)  1.43  1.43       1.4"!       1.43       1.43       1.43       1.43 

1.29       1.29  .93  1.17  1.05          .80       1.15       1.17       1.15      1.12 


WF 
max 

7600 

WF  . 
min 

900 

3 
max 

40 

min 

0 

WB 
max 

WB  . 
min 

0 

P/D 
'  max 

1.43 

P/D  ■ 
'  mm 

.01 

WF 

+  3200 

8 

+  59 

WB 

P/D 

+  .12 

30.5   34.7   1090   1920    611  (47  7  0)  (4  6  50)  1330  (3  310)  3130 
-45.3  -43.3  -10.2  -12.7  -12.7   -1.3  -17.2  -16.5   -3.7   -6.2 

11.0   19.1   2.58    246   constant  = 

(2^59)  (l7b~2)  (T .'SO) £2 '.91^ 


NOTE:   Circled  items  indicate  constraint  limit  exceeded.   The  design  code 
used  above  is  explained  as  follows : 

(1st  Digit)       -       (2nd  Digit)  (letter) 

+                       t  t 

Design  form             Power  Code  Variation  Code 

from  Table  I            1  =  70-751  from  next  page 

2  =  30-85% 
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TABLE  IV  (Cont'd! 


Variation  Code  (Also  see  Table  II] 


ABC 


Bleed 

Variable        X    X 

Constant  X    X    X    X 

State  Variables 

Scaled  X    X    X 

Unsealed         X    X    X 

Q  and  R  Weights 

Variable        X    X    X    X    X 

r.  Fixed  X 

Ship  Speed 

Constrained  XXX 

Unconstrained    X  XX 
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TABLE  V 

COMPARISON  OF  INITIAL  CONDITION  ACTUAL  OUTPUTS  WITH  THOSE 
COMPUTED  BY  THE  LINEAR  MODEL 


Output 
SMC 
T4 
QPT 
P3 

HPPT 
T5 


Obtained  from  data  generated  in  conjunction  with  Ref .  2 


Actual  Vs 

tl 

1 

ue 

Computed  Value 

32.39 

30.00 

1776 

1353 

9706 

8182 

98.53 

91.09 

3219 

-1795 

1316 

1365 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

CONCLUSIONS 

1.  Some  form  of  state  variable  scaling  is  necessary  to 
improve  the'  efficiency  of  the  numerical  calculations. 
Scaling  causes  the  A,  B,  C,  D  matrix  elements  to  be 
much  closer  to  the  same  order  of  magnitude. 

2.  If  bleed  is  to  be  utilized  as  a  control  variable, 
the  failure  of  the  mathematical  model  to  prevent  it 
from  assuming  negative  values  must  be  corrected. 

3.  Although  the  system's  non-linearity  was  in  evidence 
at  the  initial  condition  state  as  shown  in  Table  V , 
the  mathematical  model  became  a  much  more  accurate 
representation  of  the  actual  system  as  the  final 
steady  state  operating  point  was  approached. 

4 .  Optimal  Controller  design  with  the  CONSYN  program  is 
by  no  means  a  completely  automated  (or  cookbook) 
process.   Much  manual  interaction  with  options  and 
variable  parameters  must  be  employed  in  order  to 
achieve  a  feasible  design. 

5.  Input  and  input  rate  constraints  were  design  limit- 
ing.  As  long  as  these  constraints  were  not  violated, 
there  were  no  unreasonable  excursions  of  any  state  or 
output  variables. 
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B.   RECOMMENDATIONS 

1.   The  specific  reason (s)  why  the  optimization  routine 

is  so  frequently  unable  to  reduce  violated  constraints 
beyond  a  certain  point  must  be  investigated  and 
corrected  before  any  practical  controller  design  by 
this  method  can  be  achieved. 
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APPENDIX  A 
LIST  OF  CONSYN  SUBPROGRAMS 

Name  Purpose 

ANALIZ     Integral  control  design  and  analysis  (pre- 
viously named  CYCLE  in  Ref.  1) 

CONTL      Controllability  Check 

COST       Evaluation  of  Performance  index  and  deter- 
mination of  maximum  eigenvalue 

DET       Matrix  determinant 

INVERT     Matrix  inverse 

LIAP       Solution  to  bilinear  symmetric  Lyapunov 
equations 

MULT  Matrix  multiplication 

MULTV  Vector  multiplication 

MULTZ  Complex  matrix  multiplication 

OBSERV  Observability  check 

OPGAIN  Solution  to  optimal  linear  regulator 

PEAK  System  simulation 

PHIOFT  State  transition  matrix 

POLY  Coefficients  of  characteristic  polynomial 

PRINTM  Matrix  print,  double  precision 

PRINTN  Matrix  print,  single  precision 

RANK  Matrix  rank 

ROOT  Polynomial  roots 
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APPENDIX  B 


PROPULSION  PLANT  DATA 


1.  State  Variable  Representation  at  70%,  75%, 
8  0%,  and  8  5%  of  Maximum  Gas  Generator  Speed 


70%      75%      80%      35% 


States 

NG  (rpm)  6879.    7370.    7362.    3353. 

NPT  (rpm)  1058.    1313.    1742.    2300. 

V  (knots)  9.688    11.91    15.64    20.40 

Inputs 

WF  (lbm/hr)  1319.  1640.  2361.  3857. 

3  (deg)  33.29  30.53  23.37  13.20 

WB  (lbm/S)  .02  .02  .02  .02 

P/D  (••')  1.43  1.43  1.43  1.43 

Outputs 

SMC  (%) 
T4  (°R) 
QPT  (ft-lbf) 
P3  (PSIA) 
HPPT  (hp) 
T5  (°R) 


32.63 

33.21 

32.39 

27.76 

1613. 

1663. 

1776. 

2002. 

3745. 

5627. 

9706. 

17040 

54.02 

69.24 

98.53 

144  .6 

754  .4 

1407. 

3219, 

7461. 

1284  . 

1277. 

1316. 

1456. 
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APPENDIX  B  (Cont'd) 


2.  A,  B,  C,  D  MATRIX  ELEMENTS  FOR  7  5% 
OF  MAXIMUM  GAS  GENERATOR  SPEED 


The  A  Matrix 

-.2709      -.1186      0.0 
.1295      -.9364      57.71 
0.0        .0002429    -.02793 

The  B  Matrix 


.7287 

37.74 

-28.29 

0.0 

.08656 

-8 .978 

-7.861 

-526.5 

0.0 

0.0 

0.0 

.1114 

The  C  Matrix 

.01078 

-.2160 

3.171 

.02112 

.7928 

-.2289 

The  D  Matrix 


-.0003271 

0. 

.0 

.0008184 

0. 

.0 

-2.473 

0. 

.0 

.0001619 

0, 

.0 

.4526 

0. 

.0 

.01866 

0. 

.0 

-.01811 

-.5678 

2.453 

0.0 

.4905 

9.434 

10.21 

0.0 

2.120 

-219.9 

-192.5 

0.0 

.008820 

-1.405 

-1.195 

0.0 

.5300 

-54.97 

-48.13 

0.0 

.3635 

11.49 

11.85 

0.0 
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APPENDIX  B  (Cont'd) 


3.  A,  B,  C,  D  MATRIX  ELEMENTS  FOR  3  5% 
OF  MAXIMUM  GAS  GENERATOR  SPEED 


The  A  Matrix 

-1.212      -.1923      0.0 
.3466       -1.711      100.3 
0.0  .0004290    -.05131 

The  B  Matrix 


.6284 

68.15 

-26.72 

0.0 

.1194 

-14.23 

-10.79 

-1634 

0.0 

0.0 

0.0 

.3445 

The  C  Matrix 


.01097 

0.0 

0.0 

-.2715 

0.0 

0.0 

8  .489 

-5.786 

0.0 

.04678 

0.0 

0.0 

3.717 

7.101 

0.0 

-.2484 

.01736 

0.0 

The  D  Matrix 


-.007521 

-.3569 

1.229 

0.0 

.2423 

7.875 

7.273 

0.0 

2.925 

-348.4 

-264.3 

0.0 

.008168 

-1.302 

-1.349 

0.0 

1.281 

-152.6 

-115.8 

0.0 

.1855 

7.277 

6.693 

0.0 
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APPENDIX  B  (Cont'd] 


4.  SCALED  A,  B,  C,  D  MATRIX  ELEMENTS  FOR  3  5% 
MAXIMUM  GAS  GENERATOR  SPEED,  CONSTANT  BLEED 


The  Scaled  A  Matrix 

-1.212      -.2185  0.0 

.305        -1.711  .3556 

0.0         .05029  -.05181 

The  Scaled  B  Matrix 

1.915       1.412  0.0 

.3201       -.2594  -2.928 

0.0         0.0  .07237 

The  Scaled  C  Matrix 


1.163 

0  .0 

0  .0 

-.5899 

0  .0 

0.0 

.5683 

-.4402 

0.0 

.4986 

0  .0 

0  .0 

.4302 

.9341 

0  .0 

-.8712 

.06919 

0  .0 

The  Scaled  D  Matrix 


-2.43 

-.7839 

0  .0 

1.604 

.3544 

0  -0 

.5966 

-.4831 

o.o 

.2652 

-.3978 

0.0 

.4518 

-.3659 

0  .0 

1.982 

.5286 

0  .0 
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APPENDIX  B  (Cont'd) 
5.  PHYSICAL  CONSTRAINTS 


Note 


Inputs 

900  <  WF  <  7600  lbm/hr  (1) 

0  <  3  <  40°  (1) 

.01  <  P/D  <  1.43  (2) 

Input  Rates 

-3200  <  WF  <  3200  Ibm/hr/S  (3) 

-59  <  3  <  59  °/S  (3) 

-.12  <  P/D  <  .12  1/S  (4) 


NOTES:  (1)   Approx.  range  between  idle  and  full  power 
condition . 

(2)  Model  does  not  allow  for  negative  pitch. 

(3)  Obtained  from  G.E. 

(4)  Computed  from  data  provided  by  Bird  Johnson 
Company  as  follows: 

Propeller  Diameter     =   16.5  ft 
Maximum  Ahead  Pitch    =   23.5  ft 
Maximum  Reverse  Pitch  =   14.7  ft 
Minimum  Time  for  Pitch  Change  from 
Full  Reverse  to  Full  Ahead  =  2  0  sec 

p/D  =  ^W   =   23^5  _    14,7   =   12  1/s 

At     16  .  5 16  .  3 

20  S 
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APPENDIX  C 


DERIVATION  OF  EQUATIONS  USED 
FOR  SYSTEM  SIMULATION  WITH  CSMP 


The  following  matrix  equations  may  be  obtained  from 
inspection  of  the  system  block  diagram  in  Fig.  Id: 

X  =  AX  +  BU 

Y  =  CX  +  DU 

U  =  /[H(Z-EY) ]dt  -  LX 

or,  upon  differentiating, 

U  =  H(Z-EY)  -  LX 
However,  the  demand  vector  Z  is  equal  to  zero  since  all 
desired  output  deviations  are  zero  at  the  final  steady 
state  operating  point  which  has  been  defined  as  the 
origin .   Thus , 

U  =  -(HEY  +  LX) 
Recalling  that  the  matrix  E  was  defined  as 

E  *  [I  JO], 
the  above  matrix  equations  can,  therefore,  be  expanded 
into  the  following  form  for  the  3-state,  3- input,  and 
8-output  system  considered: 

*1  =  *11  Xl  +  a12  X2  +  a13  X3  +  bll  Ul  +  b12  U2  +  b13  U3 
x3  =  a31  xL  4-  a32  x2  +  a33  x3  +  b31  u±    +b32  u2  -   b33  U3 


yl  =  °11  Xl  +  C12  x2  +  C13  x3  +  dll  Ul  +  d12  u2  +  d13  U3 

Y«      =     CQ1       Xl       +     CQO     X9      +     CQ1      **      +     dQ1       U,       +     dQ9     U„      +     dQ^      U. 
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APPENDIX  C  (Cont'd) 


Ul  =  "(1?11  yl    +  hL2  72  +  h13  ^3  +  Lll  Xl  +  L12  X2  +  X13  X3} 
^3  =  _(h31  *1  +  h32  y2  +  h33  y3  +  L31  Xl  +  X32  x2  +  L33  X3} 
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APPENDIX  D 

RATIONALE  FOR  DESIGN  FORM  SELECTION 

In  Ref.  5,  the  following  outputs  were  selected  for 
feedback  in  the  case  of  the  F4  01  turbofan  je-c  engine 
controller  design: 

1.  Fan  inlet  guide  vane  angle  (FIGV) 

2.  Rear  compressor  variable  vane  angle  (RCW) 

3.  Thrust 

4.  High  turbine  inlet  temperature  (HTIT) 

Those  outputs  selected  for  feedback  in  design  form  numbe: 
2  of  this  study  were: 

1.  Propeller  pitch  to  diameter  ratio  (?/D) 

2.  Compressor  stator  guide  vane  angle  (3) 

3.  Power  turbine  torque  (QPT) 

4  .   Power  turbine  inlet  temperature  (T5) 
A  comparison  of  these  designs  follows: 

1.  The  marine  propulsion  gas  turbine  engine,  of 
course,  has  no  turbofan.   However,  the  P/D  was 
considered  analogous  to  the  FIGV  since  both 
directly  relate  to  the  loading  of  the  output 
turbine.   In  the  aircraft  case,  the  output; 
turbine  drives  the  fan;  in  the  marine  case, 
the  output  turbine  drives  the  propeller. 

2.  3  is  completely  analogous  to  RCW  since  this  is 
the  compressor  vane  control  in  both  cases. 


48 


APPENDIX  D  (Cont'd) 

The  thrust  in  the  aircraft  case  is  a  measure  of 
the  total  power  produced.   In  the  marine  case, 
either  ship  speed,  power  turbine  power,  or  power 
turbine  torque  is  analogous.   Torque  was  selected 
since  it  was  considered  the  easiest  of  the  three 
to  measure  accurately. 

HTIT  is " equivalent  to  combustor  outlet  temperature 
However,  the  high  temperature  involved  shortens 
transducer  lifetime,  and  is  variable  around  the 
perimeter  of  the  combustor.   To  eliminate  these 
difficulties,  power  turbine  inlet  temperature, 
instead,  was  chosen  for  feedback. 
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APPENDIX  E 


DERIVATION  OF  SCALED  STATE  VARIABLE  DATA 
WITH  CONSTANT  BLEED 


From  the  matrix  equation, 

X  =  AX  +  BU 

consider  one  of  the  scalar  equations 

*1    =  all  Xl  +  a12  x2  +  a13  X3  +  bll  Ul  +  b12  U2 

+  b13  u3  4-  b14  u4 

'  s       '  s 
Since  the  x    and  u    are  deviations  from  the  steady 

state  operating  point,  and  bleed  is  constant,  then  u-,  =  0 

Now  define  the  following  scaled  variables: 
x  •  u  . 

X.   = ,    U.   =  7 ,    U„  =  U„ 

x    Ax.      i    Au.      4     4 
x  x 

where  the  A  3  refer  to  the  difference  between  initial  and 

final  conditions,  and  u.    will  not  be  scaled  since  its 

initial  and  final  values  are  equal. 

Then 

AX1X1  =  all  AX1  Xl  +  a12  Ax2  X2  +  a13  "X3  X3 

+  bL1  AuL  u1  +  b12  Au2  u2  +  b14  u4 

Dividing  through  by  Ax,  results  in 

•  Ax2  Ax3  A        Au-, 

XI  =  all  Xl  +  a12  Ax^  X2  +  a13  AxY  X3  +  bll  Ax^  *1 

Au2  -     b14  - 

+  b,  0  -r un  + u, 

12  Ax-,   2    Ax-,   4 
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APPENDIX  E  (Cont'd) 

Carrying  this  procedure  out  with  the  other  state  and 
output  equations  results  in  the  scaled  matrix  equations 

X  =  AX  +  BU 
Y  =  CX  +  DU 


where 


A  = 


lll 


Ax 


1 


•21  Ax, 


Ax, 


'31  Ax. 


Ax. 

L12  Ix~ 


'22 


Ax2 
'32  Ax. 


Ax. 
■13  a3T 


Ax. 


23  Ax. 


a-, -j 

3  J 


B  = 


Au_3 
511  Ax, 


D 


All, 

21  AX 


Au 
512  Ax 


2 


1 

Au9 

*22  AxT 


14  Ax, 


!4  Ax. 


Au 
'31  AX 


1 


Au. 
'32  Ix. 


34  Ax. 
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APPENDIX    E     (Cont'd) 


C    = 


Ax, 
:11    Ay^ 

Ax-, 


'21    Ay2 

Ax, 


Ax2 
■12    Ay^ 

Ax„ 


■22    Ay2 

Ax„ 


•61    Ay6  62    Ay6 


Ax. 
'13   Ay^ 


Ax. 
''23    Ay^ 


Ax_ 
'63    W~6 


D    = 


Au-, 
lll    AyT 


Au- 


Au. 


21    Ay. 


Au, 


12    Ay 


1 


Au2 
'22    A7J 

Au0 


61    *Y6  62    ^6 


•14    Ay, 


•24    Ay, 
•64    AyJ 
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APPENDIX    F 


COMPUTER   OUTPUT    OF    OPTIMIZATION    RUN    FOR    DESIGN    6-2F 
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